Leveraging deep learning algorithms for synthetic data generation to design and analyze biological networks.

The use of synthetic data is gaining an increasingly prominent role in data and machine learning workflows to build better models and conduct analyses with greater statistical inference. In the domains of healthcare and biomedical research, synthetic data may be seen in structured and unstructured formats. Concomitant with the adoption of synthetic data, a sub-discipline of machine learning known as deep learning has taken the world by storm. At a larger scale, deep learning methods tend to outperform traditional methods in regression and classification tasks. These techniques are also used in generative modeling and are thus prime candidates for generating synthetic data in both structured and unstructured formats. Here, we emphasize the generation of synthetic data in healthcare and biomedical research using deep learning methods for unstructured data formats such as text and images. Deep learning methods leverage the neural network algorithm, and in the context of generative modeling, several neural network architectures can create new synthetic data for a problem at hand including, but not limited to, recurrent neural networks (RNNs), variational autoencoders (VAEs), and generative adversarial networks (GANs). To better understand these methods, we will look at specific case studies such as generating realistic clinical notes of a patient, the generation of synthetic DNA sequences, as well as to enrich experimental data collected during the study of heterotypic cultures of cancer cells.

Intrinsically Disordered Proteins: Critical Components of the Wetware.

Despite the wealth of knowledge gained about intrinsically disordered proteins (IDPs) since their discovery, there are several aspects that remain unexplored and, hence, poorly understood. A living cell is a complex adaptive system that can be described as a wetware─a metaphor used to describe the cell as a computer comprising both hardware and software and attuned to logic gates─capable of "making" decisions. In this focused Review, we discuss how IDPs, as critical components of the wetware, influence cell-fate decisions by wiring protein interaction networks to keep them minimally frustrated. Because IDPs lie between order and chaos, we explore the possibility that they can be modeled as attractors. Further, we discuss how the conformational dynamics of IDPs manifests itself as conformational noise, which can potentially amplify transcriptional noise to stochastically switch cellular phenotypes. Finally, we explore the potential role of IDPs in prebiotic evolution, in forming proteinaceous membrane-less organelles, in the origin of multicellularity, and in protein conformation-based transgenerational inheritance of acquired characteristics. Together, these ideas provide a new conceptual framework to discern how IDPs may perform critical biological functions despite their lack of structure.

Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma.

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and "promiscuous" interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discerning the energy landscape representation of the weakly funneled IDPs in terms of reaction coordinates is challenging. To understand conformational dynamics in real time and decipher how IDPs recognize multiple binding partners with high specificity, several sophisticated knowledge-based and physics-based in silico sampling techniques have been developed. Here, using specific examples, we highlight recent advances in energy landscape visualization and molecular dynamics simulations to discern conformational dynamics and discuss how the conformational preferences of IDPs modulate their function, especially in phenotypic switching. Finally, we discuss recent progress in identifying small molecules targeting IDPs underscoring the potential therapeutic value of IDPs. Understanding structure and function of IDPs can not only provide new insight on cellular decision making but may also help to refine and extend Anfinsen's structure/function paradigm.

Group Behavior and Emergence of Cancer Drug Resistance.

Drug resistance is a major impediment in cancer. Although it is generally thought that acquired drug resistance is due to genetic mutations, emerging evidence indicates that nongenetic mechanisms also play an important role. Resistance emerges through a complex interplay of clonal groups within a heterogeneous tumor and the surrounding microenvironment. Traits such as phenotypic plasticity, intercellular communication, and adaptive stress response, act in concert to ensure survival of intermediate reversible phenotypes, until permanent, resistant clones can emerge. Understanding the role of group behavior, and the underlying nongenetic mechanisms, can lead to more efficacious treatment designs and minimize or delay emergence of resistance.

Dynamic Phenotypic Switching and Group Behavior Help Non-Small Cell Lung Cancer Cells Evade Chemotherapy.

Drug resistance, a major challenge in cancer therapy, is typically attributed to mutations and genetic heterogeneity. Emerging evidence suggests that dynamic cellular interactions and group behavior also contribute to drug resistance. However, the underlying mechanisms remain poorly understood. Here, we present a new mathematical approach with game theoretical underpinnings that we developed to model real-time growth data of non-small cell lung cancer (NSCLC) cells and discern patterns in response to treatment with cisplatin. We show that the cisplatin-sensitive and cisplatin-tolerant NSCLC cells, when co-cultured in the absence or presence of the drug, display dynamic group behavior strategies. Tolerant cells exhibit a 'persister-like' behavior and are attenuated by sensitive cells; they also appear to 'educate' sensitive cells to evade chemotherapy. Further, tolerant cells can switch phenotypes to become sensitive, especially at low cisplatin concentrations. Finally, switching treatment from continuous to an intermittent regimen can attenuate the emergence of tolerant cells, suggesting that intermittent chemotherapy may improve outcomes in lung cancer.

Protein conformational dynamics and phenotypic switching.

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure but exist as conformational ensembles. Because of their structural plasticity, they can interact with multiple partners. The protein interactions between IDPs and their partners form scale-free protein interaction networks (PINs) that facilitate information flow in the cell. Because of their plasticity, IDPs typically occupy hub positions in cellular PINs. Furthermore, their conformational dynamics and propensity for post-translational modifications contribute to "conformational" noise which is distinct from the well-recognized transcriptional noise. Therefore, upregulation of IDPs in response to a specific input, such as stress, contributes to increased noise and, hence, an increase in stochastic, "promiscuous" interactions. These interactions lead to activation of latent pathways or can induce "rewiring" of the PIN to yield an optimal output underscoring the critical role of IDPs in regulating information flow. We have used PAGE4, a highly intrinsically disordered stress-response protein as a paradigm. Employing a variety of experimental and computational techniques, we have elucidated the role of PAGE4 in phenotypic switching of prostate cancer cells at a systems level. These cumulative studies over the past decade provide a conceptual framework to better understand how IDP conformational dynamics and conformational noise might facilitate cellular decision-making.

Role of immunotherapy and co-mutations on KRAS-mutant non-small cell lung cancer survival.

BACKGROUND: KRAS mutations reported in non-small cell lung cancer (NSCLC) represent a significant percentage of patients diagnosed with NSCLC. However, there still remains no therapeutic option designed to target KRAS. In an era with immunotherapy as a dominant treatment option in metastatic NSCLC, the role of immunotherapy in. KRAS: mutated patients is not clear. METHODS: Eligible patients diagnosed with NSCLC and found to have a KRAS mutation were identified in an institutional lung cancer database. Demographic, clinical, and molecular data was collected and analyzed. RESULTS: A total of 60 patients were identified for this retrospective analysis. Majority of patients were Caucasian (73%), diagnosed with stage IV (70%) adenocarcinoma (87%), and had a KRAS codon 12 mutation (78%). Twenty percent of patients were treated with immunotherapy. Median overall survival was 28 months in the cohort and patients who received immunotherapy were found to have better survival versus those who did not (33 vs. 22 months, P=0.31). Furthermore, there was an association between high survival and patients who received immunotherapy (P=0.007). CONCLUSIONS: Patients with KRAS mutations have a unique co-mutation phenotype that requires further investigation. Immunotherapy seems to be an effective choice of treatment for KRAS positive patients in any treatment-line setting and yields better outcomes than conventional chemotherapy. The relationship between immunotherapy and KRAS mutations requires further studies to confirm survival advantage.

Machine Learning for Prioritization of Thermostabilizing Mutations for G-Protein Coupled Receptors.

Although the three-dimensional structures of G-protein coupled receptors (GPCRs), the largest superfamily of drug targets, have enabled structure-based drug design, there are no structures available for 87% of GPCRs. This is due to the stiff challenge in purifying the inherently flexible GPCRs. Identifying thermostabilized mutant GPCRs via systematic alanine scanning mutations has been a successful strategy in stabilizing GPCRs, but it remains a daunting task for each GPCR. We developed a computational method that combines sequence-, structure-, and dynamics-based molecular properties of GPCRs that recapitulate GPCR stability, with four different machine learning methods to predict thermostable mutations ahead of experiments. This method has been trained on thermostability data for 1231 mutants, the largest publicly available data set. A blind prediction for thermostable mutations of the complement factor C5a receptor 1 retrieved 36% of the thermostable mutants in the top 50 prioritized mutants compared to 3% in the first 50 attempts using systematic alanine scanning.

Monitoring and Determining Mitochondrial Network Parameters in Live Lung Cancer Cells.

Mitochondria are dynamic organelles that constantly fuse and divide, forming dynamic tubular networks. Abnormalities in mitochondrial dynamics and morphology are linked to diverse pathological states, including cancer. Thus, alterations in mitochondrial parameters could indicate early events of disease manifestation or progression. However, finding reliable and quantitative tools for monitoring mitochondria and determining the network parameters, particularly in live cells, has proven challenging. Here, we present a 2D confocal imaging-based approach that combines automatic mitochondrial morphology and dynamics analysis with fractal analysis in live small cell lung cancer (SCLC) cells. We chose SCLC cells as a test case since they typically have very little cytoplasm, but an abundance of smaller mitochondria compared to many of the commonly used cell types. The 2D confocal images provide a robust approach to quantitatively measure mitochondrial dynamics and morphology in live cells. Furthermore, we performed 3D reconstruction of electron microscopic images and show that the 3D reconstruction of the electron microscopic images complements this approach to yield better resolution. The data also suggest that the parameters of mitochondrial dynamics and fractal dimensions are sensitive indicators of cellular response to subtle perturbations, and hence, may serve as potential markers of drug response in lung cancer.

Modeling small cell lung cancer (SCLC) biology through deterministic and stochastic mathematical models.

Mathematical cancer models are immensely powerful tools that are based in part on the fractal nature of biological structures, such as the geometry of the lung. Cancers of the lung provide an opportune model to develop and apply algorithms that capture changes and disease phenotypes. We reviewed mathematical models that have been developed for biological sciences and applied them in the context of small cell lung cancer (SCLC) growth, mutational heterogeneity, and mechanisms of metastasis. The ultimate goal is to develop the stochastic and deterministic nature of this disease, to link this comprehensive set of tools back to its fractalness and to provide a platform for accurate biomarker development. These techniques may be particularly useful in the context of drug development research, such as combination with existing omics approaches. The integration of these tools will be important to further understand the biology of SCLC and ultimately develop novel therapeutics.

A modified Stokes-Einstein equation for Aß aggregation.

BACKGROUND: In all amyloid diseases, protein aggregates have been implicated fully or partly, in the etiology of the disease. Due to their significance in human pathologies, there have been unprecedented efforts towards physiochemical understanding of aggregation and amyloid formation over the last two decades. An important relation from which hydrodynamic radii of the aggregate is routinely measured is the classic Stokes-Einstein equation. Here, we report a modification in the classical Stokes-Einstein equation using a mixture theory approach, in order to accommodate the changes in viscosity of the solvent due to the changes in solute size and shape, to implement a more realistic model for Aß aggregation involved in Alzheimer's disease. Specifically, we have focused on validating this model in protofibrill lateral association reactions along the aggregation pathway, which has been experimentally well characterized. RESULTS: The modified Stokes-Einstein equation incorporates an effective viscosity for the mixture consisting of the macromolecules and solvent where the lateral association reaction occurs. This effective viscosity is modeled as a function of the volume fractions of the different species of molecules. The novelty of our model is that in addition to the volume fractions, it incorporates previously published reports on the dimensions of the protofibrils and their aggregates to formulate a more appropriate shape rather than mere spheres. The net result is that the diffusion coefficient which is inversely proportional to the viscosity of the system is now dependent on the concentration of the different molecules as well as their proper shapes. Comparison with experiments for variations in diffusion coefficients over time reveals very similar trends. CONCLUSIONS: We argue that the standard Stokes-Einstein's equation is insufficient to understand the temporal variations in diffusion when trying to understand the aggregation behavior of Aß42 proteins. Our modifications also involve inclusion of improved shape factors of molecules and more appropriate viscosities. The modification we are reporting is not only useful in Aß aggregation but also will be important for accurate measurements in all protein aggregation systems.

Stability of two cluster solutions in pulse coupled networks of neural oscillators.

Phase response curves (PRCs) have been widely used to study synchronization in neural circuits comprised of pacemaking neurons. They describe how the timing of the next spike in a given spontaneously firing neuron is affected by the phase at which an input from another neuron is received. Here we study two reciprocally coupled clusters of pulse coupled oscillatory neurons. The neurons within each cluster are presumed to be identical and identically pulse coupled, but not necessarily identical to those in the other cluster. We investigate a two cluster solution in which all oscillators are synchronized within each cluster, but in which the two clusters are phase locked at nonzero phase with each other. Intuitively, one might expect this solution to be stable only when synchrony within each isolated cluster is stable, but this is not the case. We prove rigorously the stability of the two cluster solution and show how reciprocal coupling can stabilize synchrony within clusters that cannot synchronize in isolation. These stability results for the two cluster solution suggest a mechanism by which reciprocal coupling between brain regions can induce local synchronization via the network feedback loop.

Synaptic and intrinsic determinants of the phase resetting curve for weak coupling.

A phase resetting curve (PRC) keeps track of the extent to which a perturbation at a given phase advances or delays the next spike, and can be used to predict phase locking in networks of oscillators. The PRC can be estimated by convolving the waveform of the perturbation with the infinitesimal PRC (iPRC) under the assumption of weak coupling. The iPRC is often defined with respect to an infinitesimal current as z(i)(ϕ), where ϕ is phase, but can also be defined with respect to an infinitesimal conductance change as z(g)(ϕ). In this paper, we first show that the two approaches are equivalent. Coupling waveforms corresponding to synapses with different time courses sample z(g)(ϕ) in predictably different ways. We show that for oscillators with Type I excitability, an anomalous region in z(g)(ϕ) with opposite sign to that seen otherwise is often observed during an action potential. If the duration of the synaptic perturbation is such that it effectively samples this region, PRCs with both advances and delays can be observed despite Type I excitability. We also show that changing the duration of a perturbation so that it preferentially samples regions of stable or unstable slopes in z(g)(ϕ) can stabilize or destabilize synchrony in a network with the corresponding dynamics.

Pulse coupled oscillators and the phase resetting curve.

Limit cycle oscillators that are coupled in a pulsatile manner are referred to as pulse coupled oscillators. In these oscillators, the interactions take the form of brief pulses such that the effect of one input dies out before the next is received. A phase resetting curve (PRC) keeps track of how much an input advances or delays the next spike in an oscillatory neuron depending upon where in the cycle the input is applied. PRCs can be used to predict phase locking in networks of pulse coupled oscillators. In some studies of pulse coupled oscillators, a specific form is assumed for the interactions between oscillators, but a more general approach is to formulate the problem assuming a PRC that is generated using a perturbation that approximates the input received in the real biological network. In general, this approach requires that circuit architecture and a specific firing pattern be assumed. This allows the construction of discrete maps from one event to the next. The fixed points of these maps correspond to periodic firing modes and are easier to locate and analyze for stability compared to locating and analyzing periodic modes in the original network directly. Alternatively, maps based on the PRC have been constructed that do not presuppose a firing order. Specific circuits that have been analyzed under the assumption of pulsatile coupling include one to one lockings in a periodically forced oscillator or an oscillator forced at a fixed delay after a threshold event, two bidirectionally coupled oscillators with and without delays, a unidirectional N-ring of oscillators, and N all-to-all networks.

Phase-resetting curves determine synchronization, phase locking, and clustering in networks of neural oscillators.

Networks of model neurons were constructed and their activity was predicted using an iterated map based solely on the phase-resetting curves (PRCs). The predictions were quite accurate provided that the resetting to simultaneous inputs was calculated using the sum of the simultaneously active conductances, obviating the need for weak coupling assumptions. Fully synchronous activity was observed only when the slope of the PRC at a phase of zero, corresponding to spike initiation, was positive. A novel stability criterion was developed and tested for all-to-all networks of identical, identically connected neurons. When the PRC generated using N-1 simultaneously active inputs becomes too steep, the fully synchronous mode loses stability in a network of N model neurons. Therefore, the stability of synchrony can be lost by increasing the slope of this PRC either by increasing the network size or the strength of the individual synapses. Existence and stability criteria were also developed and tested for the splay mode in which neurons fire sequentially. Finally, N/M synchronous subclusters of M neurons were predicted using the intersection of parameters that supported both between-cluster splay and within-cluster synchrony. Surprisingly, the splay mode between clusters could enforce synchrony on subclusters that were incapable of synchronizing themselves. These results can be used to gain insights into the activity of networks of biological neurons whose PRCs can be measured.

Backbone structure of the amantadine-blocked trans-membrane domain M2 proton channel from Influenza A virus.

Amantadine is known to block the M2 proton channel of the Influenza A virus. Here, we present a structure of the M2 trans-membrane domain blocked with amantadine, built using orientational constraints obtained from solid-state NMR polarization-inversion-spin-exchange-at-the-magic-angle experiments. The data indicates a kink in the monomer between two helical fragments having 20 degrees and 31 degrees tilt angles with respect to the membrane normal. This monomer structure is then used to construct a plausible model of the tetrameric amantadine-blocked M2 trans-membrane channel. The influence of amantadine binding through comparative cross polarization magic-angle spinning spectra was also observed. In addition, spectra are shown of the amantadine-resistant mutant, S31N, in the presence and absence of amantadine.